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ABSTRACT 


A method for remote sensing of cloud parameters by using infrared 
sounder data has been developed on the basis of the parameterized infrared 
transfer equation applicable to cloudy atmospheres. The method is utilized 
for the retrieval of the cloud height, amount, and emlssivity in 11 ym 
region. Numerical analyses and retrieval experiments have been carried out 
by utilizing the synthetic sounder data for the theoretical study. The 
sensitivity of the numerical procedures to the measurement and instrument 
errors are also examined. The retrieved results are physically discussed 
and numerically compared with the model atmospheres. Comparisions reveal 
that the recovered cloud parameters agree reasonably well with the 
pre-assumed values. However, for cases when relatively thin clouds and/or 
small cloud fractional cover within a field of view are present, the 
recovered cloud parameters show considerable fluctuations . 

Experiments on the proposed algorithm are carried out utilizing High 
Resolution Infrared Sounder (HIRS/2) data of NOAA 6 and TIROS-N. Results 
of experiments show reasonably good comparisons with the surface reports 
and GOES satellite images. 


/ 
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1. Introduction 

Satellite measurements can provide global coverage of cloud properties 
such as cloud top height and temperature, fractional cloud cover and cloud 
emisslvity. Those data are essential for the inclusion of cloud 
parameterizations in realistic models of the earth's climate (Schneider, 

1972; Cess, 1974; Stephens and Webster, 1979). Recently, the World Climate 
Research Programme (WCRP) plan recognized the need to develop a uniform 
global cloud climatology as part of a broad research program on climate 
processes. The International Satellite Cloud Climatology Project (ISCCP) 
thus has been approved as the first project of WCRP. The basic objective 
of the ISCCP Is to collect and analyze satellite radiance data to infer the 
global distribution of cloud radiative properties in order to improve the 
modeling of cloud effects on climates. 

The purpose of this research is to explore the algorithm retrieving 
the cloud properties by utilizing the available infrared sounder data from 
polar-orbiting satellites. Although the geostationary meteorological 
satellite measurements are used as the primary operational data for the 
ISCCP, the polar-orbiting satellites are essential to the project providing 
1) coverage of the high latitude regions not viewed by the geostationary 
satellites; 2) a basis for normalization of the radiances observed by the 
different geostationary satellites; 3) global coverage that may help 
mitigate the possible loss of one or more geostationary satellites; 

4) multispectral measurements for discriminating cloud properties not 
derivable from the primary two-channel geostationary data (Schiffer and 
Rossow, 1983). In addition to the HIRS data on polar-orbiting 
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satellites, the VISSR Atmosperic Sounder (VAS) data of GOES satellites 
(Smith et al., 1981) can also be applied to the current proposed algorithm 
in future. 

In this study, two channels of 15 um CC,. assorption band and one 
channel of 11 pm window region are used for the retrieval of cloud top 
height, fractional cloud cover and emissivity of 11 ym region. The 
application of sounder data to cloud parameters retrieval has been 
investigated previously. McCleese and Wilson (1976), Smith and Platt 
(1978) and Wielicki and Coakley (1981) have demonstrated that cloud heights 
and effective cloud cover can be recovered under a wide variety of 
atmospheric conditions using different regions in the near and far infrared 
parts of the terrestrial spectrum. Also, recently, Chahine (1982) applied 
15 pm infrared radiance data for day and night mapping of the global 
distributions of the horizontal cloud covers and the corresponding cloud 
top pressure levels. His results showed that the effective cloud cover 
derived from 15 pm data is less than that obtained from visible data. The 
discrepancy is obviously due to the existence of nonblack clouds in the 
atmosphere. It Is therefore important for us to determine the cloud 
emissivity and "physical cloud cover" separately in the retrieval. Only 

i 

such cloud cover can be compared with values obtained from nephanalysis of 
visible cloud images . 

In the following sections, we present the basic formulation of the 
upwelling radiance in cloudy atmospheres and the technique for the cloud 
parameters retrieval . Numerical experiments and a number of case studies 
utilizing the HIRS/2 data for the retrievals are then described and 
physically discussed. 



2. Cloud Parameters Retrieval 


Consider an atmosphere containing fractional cloud cover A within a 
satellite instantaneous field-of-view (ITOV). The emitted radiation 
Intercepted by the satellite radiometer can be written in the equation 
calculating the upwelling radiance at top of the atmosphere, 

1(^,9,-) - AI c (v i ,9,«)+(l-A)I(u i> 9,-), (1) 

where is the dear-column mean spectral upwelling radiance 

measured in a channel whose mean effective wavelength is and the local 
zenith angle of the observation is 9. If we consider only the thermal 
emission in the atmosphere at longwave infrared frequencies, 1(^,9,®) can 
be described by the solution of the radiative transfer equation under the 
condition of thermodynamic equilibrium, 

00 

I(v 1 ,9,«) - B(w 1 ,T s )T(w 1 ,9;-,0)+/ B[v i ,T(Z)IK(v i ,9;-, Z)dZ, (2) 

o 

where B is the spectral Planck function associated with temperature T, T is 
the spectral transmittance for the relevant gas, is surface temperature, 
and K( 9;®,Z)=dT( v^, 9;®,2)/dz is the spectral weighting function. In 
Eq. (1), I ( y ^,9,®) is the upwelling radiance emitted from cloudy area, and 
can be derived in terms of the cloud and atmospheric radiative properties. 
This radiance is caused by I) the transmission of the upwelling radiance 
from the cloud layer, 2) the emission from the cloud top, and 3) the 
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emission and absorpcion contribution of the gases above the cloud top. 

Their mathematical expressions can be written in the following forms, 

I (^,8,-) = {l(v i ,8,Z t )T c (v i ,8) 

+B(v 1> T(Z)l[l-T c (v i ,8)]}T(v i ,0;« Z t ) 

+/“b[v T(Z)K(v ,8;-,Z)dZ, (3) 

2 t 

where Z £ is cloud top height, and the cloud transmissivity which 

is function of the cloud thickness and optical properties . The upwelling 
radiance from the atmosphere below the - cloud layer I(\^,8,Z c ) can be 
expressed in terms of the radiative transfer equation 

I(V8,Z e ) * B < v i .T s )T(v i ,8;Z t ,0) 

, Z t 

+/B[v T(Z)]K(v,8;Z,Z)dZ. (4) 

o * 1 

In writing the preceding equations, the surface infrared emissivity is 
assumed to be unity, and the cloud reflectivity is considered to be 
negligible (e*g., Paltridge and Platt, 1976). The cloud radiative 

o 

properties are described by the cloud transmissivities which depend on the 
temperature structure and the total water path in the vertical column of 
the cloud. All quantities given in the above equations are integrated over 
wavenumber, and weighted by the spectral response of the instrument. 
However, since B varies slowly with v, while T varies rapidly and without 
correlation to B within the narrow spectral channels of the HIB.S/2, it 
suffices to perform the spectral integrations of B and T independently and 
treat the results as if they were monochromatic values for the effective 
wavenumber (Yeh and Liou, 1983; Smith, 1983). We thus may change the 
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reference layer of clear column transmitcances and weighting functions from 
cloud layer to the top of the atmophere in Eq . (4), 

T( ^ , 9 ; Z ) - T(v i> 0;Z t> Z)T(v i ,9;-,Z t ) (5) 

K^.Qj-.Z) - K(v i ,8;Z t ,Z)K(v 1 ,8;-,Z). 

For simplicity, we eliminate 9 in the equations, replace by i to indicate 
the channel number, and define T( v 8;°», z)=>T^(z) and K( , 8; ®, z)=«K^(z) . On 
substituting Eqs. (2), (3), (4) and (5) into (1), we obtain parameterized 
equation for the upwelling radiance emergent from a partially, single- layer 
cloudy atmosphere in the form 

- I t (-)-A(l-Tj)I^(Z t ), (6) 


where 


I i CZ t ) * B i (T s )T 1 (0)+/ o t B ;L [T(Z)]K i (z)dZ-B i (Z t; )T i (Z t; ). (7) 

In Eq. (6), two terms on the right-hand side represent, respectively, the 
dear-column radiance and the correction of the first term due to the cloud 
contamination within the IFOV. For mul tipi e-layer cloudy atmospheres, 

Eq. (6) will be used as if the lower clouds are attached to the top-layer 
cloud, and their combination can be viewed as one thick cloud. In the 

/•» 'V* 

absence of the cloud, A“0 and TJ>1, the Eq . (6) reduces to 1^ ( <*»)™I^ ( «) . 
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a. Cloud cop height determination 

Let the difference between the upwelling radiances for clear and 
cloudy atmospheres be 

- liO-V (8) 


To a good approximation, Che cloud transmissivities may be considered to be 
the same for the adjacent channels within HIE.S/2 15 urn Interval. That is, 
where i and j are two adjacent channels in the 15 um CO band. We define 
the ratio to the radiance differences for 1 and j channels as 


AG, 


H(Z C ) - 


AG. 

J 


w 

w 


(9) 


Because of the ratio, the H function is independent of the cloud opacity, 
and depends only on the weighting function of the channels and the radiance 
profile from the surface to the cloud top. Since the weighting function 
varies slightly with changing temperature and water vapor profiles , the 
approximate relation between the H value and Z £ may be constructed through 
Eq . (9) with a given radiance profile B^(Z). By evaluating the radiance 
profile from radiosonde observations or climatological profiles, we may 
determine Z £ from the value of AG.^ and AG^, which are obtained from Eq . (8) 
when observed radiance 1^ and the clear column radiance I^(®) are used. 
I^(«) may be computed based on the available radiosonde profiles. It 
should be noted that degree of the polynomial relating Z and H value may 
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affect the accuracy of the results. More detailed discussion will be given 
in the section on numerical experiments. 

The method of using adjacent channels within 15 urn region for cloud 
height inference has been discussed previously (e.g., Smith and Platt, 

1976; Wielicki and Coakley, 1981; Yeh et al., 1984). No further discussion 
about the method will be repeated here. 

(b) Cloud amount and emlsslvity determination 

The radiance of HIRS/2 channel 7 (at 13.4 inn) and channel 8 (at 11 pm) 
are mostly contributed from the gaseous emission in the lower atmosphere 
and the surface for the clear atmosphere in view of the distribution of 
weighting functions (Fig. 1). Both channels are sensitive to the presence 
of clouds at any altitude. For nonblack clouds, the cloud radiative 
properties at these two frequencies are distinct, but the relationship of 
their radiative properties can be used in determining the cloud parameters. 

Assume that the cloud transmissivity may be approximately expressed as 
an exponential function of the cloud optical depth T which is the product 
of the volume absorption coefficient 8^ and the geometric thickness 

T^(AZ) » exp(-T i ) (10) 

* exp(-8 i AZ) . 

One classical problem of computing the radiative transfer within clouds is 
defining the medium in the sense that no operational model is likely to 
predict the drop size distribution as a function of height within the 
cloud. The predictable parameters will probably be restricted to either 
liquid and/or ice water content per volume w or the total water content W 
in the vertical column of the cloud. The grossest approximation is simply 
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co assign one volume absorption coefficienc Co each cloud type, 
irrespective of the variations in cloud water content and drop size 
distributions. The volume absorption coefficient then represents a mean 
radiative property to each cloud type. 

If possible it is obviously preferable to relate the volume absorption 
coefficient as directly as possible to water content w. The ideal case is 
of a cloud which has a unique drop size distribution, since then 3^ is 
simply proportional to the total number of drops per unit volume and 
therefore to w. Paltridge and Platt (1976) found that some fairly simple 
relation does exist between 3^ and w. They showed that the value of dS^/dw 
for 11 urn is nearly a constant, and therefore 

S i - w (11) 

The parameter ^(“dS^/dw) is In fact the mass absorption coefficient by the 
above definition. Such a relation Is valid approximately for the entire 
window region of the spectrum from 8 to 14 pm. 

Considering two narrow spectra centering at 11 and 13.4 pm, the 
experiments by Feddes and Lious (1978) show that r=>K- , . /*., is 

VQS il Uni 

approximately equal to 1.1 for both ice and water clouds. If we insert Eq . 
(10) to Eq . (6), and apply the consequence to the two channels of 11 and 
13.4 pm, we may obtain the following equation 

px-qx r -(p-q) - 0, 


( 12 ) 
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where 


P - (Vvvv 

q - (I 8 -I 8 )I^(Z c ) (13) 

x - Tg. 

The subscripts 7 and 8 in the equations denote the channel number of the 
HIRS/2. It is readily seen that one may apply iterative method solving x 
in Eq . (12). Solution of x is the cloud transmissivity at 11 van, and is 
always within (0,1). The cloud eoissivity in 11 urn region can then be 
obtained by the following approximation 

e «* 1-x. (14) 

And the fractional cloud cover is 

A - tf 8 -I 8 )/eI 8 (Z c ). (15) 

One should know that solving x in Eq.(12) is not an easy task, because the 
noises from the measured radiances affecting the coefficients p and q in 
the equation. The difficulties concerning numerical calculations will be 


discussed in the next section. 



Page 1 3 


3. Numerical experiments and theoretical analyses 
a. Numerical procedures 

In this section, we construct a hypothetical cloudy atmosphere and 
perform numerical calculations of atmospheric parameters in order to test 
the validity of the foregoing derivations. The error sensitivity of the 
numerical procedures will also be examined. The measurement errors are 
introduced into theoretical calculations, and results of the errors 
generated in the sequential computations are analyzed. 

In the error analysis, we use the radiosonde profile of Charleston, 
S.C., at 1200 (all times GMT), 24 April 1980 as the model atmosphere filled 
with various types of clouds in an IFOV. The "measured" radiances are 
simulated for a range of 10 cloud amounts (0.1 to 1.0), 3 cloud top heights 
(300 mb to 700 mb) and 6 cloud optical thickness in 11 ym region (0.5 to 
3.0). In calculating the measured radiances, we introduce a Gaussian 
random noise distribution to the temperature and water vapor profiles. 

These random noises, which are to be added to 40 values of temperature and 
water vapor concentration, are determined by an inverse Gaussian 
probability distribution function. The standard deviations used for the 
temperature and water vapor profiles are l°k and 20Z, respectively. 

After this is done, we then impose some artificial errors on to the 

"measured" values to simulate instrumental noises. The noise is taken to 

- 2-1 

be Gaussian with zero mean and a standard deviation of 0.22 mWm sr cm for 
the 15 ym channels according to the noise levels specified for the HIRS/2 
Instrument on TIB.0S-N (Schwalb, 1978). For each cloud amount, optical 
depth, and cloud top height combination, 200 samples are then 
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constructed in order to investigate the statistical properties of the 
retrieved quantities. 

Before we display the results of the calculations, we shall discuss 
the difficulties of solving Eq . (12) first. Although the solution of x has 
to be within the domain of (0,1), there is a solution X s * 1 (clear 
atmosphere) always satisfies Eq. (12). While the given measurement errors, 
which affect the values p and q, cause the failure of finding a solution 
within (0,1), 3P*1 will be automatically obtained as the solution. However, 
this solution may not be acceptable for many cases. Especially when the 
cloud optical depth is large, and true solution is close to the left margin 
of the domain (0,1), one may easily fall to find a solution within (0,1). 

To minimize such mistakes, some criteria must be imposed to reexamine the 
solution if x is found to be 1. 

To illustrate the statements given above, the cloud-free (x-1) 
percentages for the various cloud amount and optical depth combinations are 
plotted and shown in Fig. 2. Without imposing any criteria, the cloud-free 
percentage obtained by solving Eq. (12) is as expected decreasing with 
increasing optical depth from 0.5 to 1.5. This is due to the fact that the 
larger cloud signatures are sensed from thicker clouds by the satellite 
radiometer than those from thinner clouds. Therefore, less probability of 
cloud-free atmosphere is obtained for thicker cloud. However, the 
cloud-free percentage rises when the optical depth is greater than 1.5. 

This is because the cloud transmissivity becomes smaller (more opaque) for 
the thicker clouds, and the same amount of measurement error may cause 
higher opportunity of failing to find a solution within (0,1). Since x=l 
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automatically becomes Che solution for such a no-solution situation, 
cloud-free percentage thus increases incorrectly. 

By examining Eq . (12), we find that for the smaller value of the 
constant term p-q, the influence from measurement errors on the solution x 
becomes greater proportionally. Recalling Eq. (13), we may have 

Apq - p-q » I^(Z t )Ig(2 c )[T^(Z t )-^(Z t )]A. (16) 

It is readily seen that Apq is linearly dependent on the cloud 
transmissivities and cloud amount. In Fig. 3, it shows that Apq generally 
decreases with increase of cloud optical depth x if t>l, but increases with 
increase of t if x<l. When a cloud is optically thick enough (or a black 
cloud), Apq should be close to 0, irrespective of the change of cloud 
amount. For clouds with the given cloud transmissivities, Apq is 0 when 
A*0, and is a maximum value when A»1 (overcast). Actually, Apq can have 
negative values when the measurement errors are involved in the cases of 
small cloud amount or optically thick cloud. (For optical very thin 
clouds, Apq may also be negative, but the solution of x should be close to 
1 anyway; no correction will be attempted in this situation.) 

If we evaluate the difference between the measured radiances of 11 and 
13.4 um channels, we find that 

“ ■ V*7 ’ 

+ pq/ig(z t ). 


( 17 ) 
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Results of Ai' are also plotted on Fig. 3. It shows that aI* has the maximum 
value when no cloud is present within IFOV, but declines with increase of 
cloud amount. In Eq. (17), when the optical depth increases, the second 

'X. 

term on the right-hand side becomes more important, and AI becomes smaller 

% 

or even a negative value. From Fig. 3, we find that the tendencies of AI 

'X, 

and Apq are oposite with variation in cloud amount. Since the value of AI 
decreases with increase of cloud amount, the differentiation of clear or 
cloudy atmospheres by comparing AI* with a criterion become feasible. Let 
us define 

R Q - <I a -I 7 )-Apq/Ig(Z t ), (18) 


then 


AI > R q , clear; (19) 


otherwise, cloudy. Note that the dimension of Apq/I„(Z ) is the same as 

o t 

% 

AI, and Apq contains the information of cloud amount and cloud optical 
depth (or emissivity) as Eq . (16) implies. In the program, Apq is obtained 
from the direct result of Eq . (13). 

To display the response of values of Apq and AI to the change of 

'X. - - 

radiosonde profile, we also plot Apq and AI with variations in cloud amount 
and optical depth for tropical and midlatitude winter climatological 

profiles (Figs. 4 and 5). The basic patterns remain the same with only 

% 

difference of the stretches of Apq and AI. These values obtained by using 
midlatitude winter climatological profile apparently are within narrower 
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Ftg. 3. Cloud amount versus aa and at (Charleston, SC) Pig. 4. Cloud Amount versus Aoo and at (Tropical climate profile) . 
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region due to colder profile and Che less distinctive upvelllng radiances 
profiles from channels 7 and 8. In contrast, the values using the tropical 
climatological profile show wider stretches of Apq and AI due to warmer 
radiosonde profile and more distinctive upwelling radiances in two 
channels . 

After reevaluating the solution of x-1 to determine if it is actually 
clear or cloudy atmosphere, it is desirable to further retrieve the cloud 
parameters based on the following procedures. It should be known that the 
results from the following discrete differentiation method is less 
sensitive to the measurement errors, although may not be as accurate as the 
solutions directly solved from Eq . (12). Again, the value of Apq can be 
used as a guideline for the criteria- If the values of Apq become 
negative, we shall expect the solutions to be large cloud optical depth (or 
emissivity) and/or small cloud amount- Since emissivity is always within 
(0,1), we decide to estimate emissivity in stead of optical depth, which, 
theoretically speaking, is between 0 and infinity. 

For the cases of Apq<0, let us assume 

R « * [<VV +4 W/VV1- S »[VV-VV]- 


m- 1, • -- ,8, 


( 20 ) 


where S^*0.Q3, S 2 =O.05, S 3 -0.09, S 4 =0.15 and S^-0.23, then 




AI>R 1 , 

r 1 >ai>r 2 , 

r 2 >ai>r 3 , 

r 3 >ai>r 4 , 

r 4 >ai>r 5 , 


e-0.20 
e- 0 . 40 
£-0 .60 
e*0.80 


( 21 ) 


e-0.95 
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Otherwise, e»1.00. 

Should Apq>0 occur, we may assume 




n**6, . . . , 10, 


( 22 ) 


where S^aO.075, S 7 *0.15, Sg=*0.20, S g =0.275 and S 1Q =0.35, then 

AI>R 6 , A-0.10; 

R 6 >AI>R 7 , A=0 . 25 ; 

R ? >AI>R 8 , A-0.40; 

R 8 >AI>R g , A-0.55; 

R 9 >/tt>R 10 , A-0.70; 


(23) 


Otherwise, A=<3.90. After solving either cloud emissivity or cloud amount 
from the above procedures, we then may go back to Eq . (6) to solve cloud 
amount or emissivity. 

It should be explained the reasons that we find cloud emissivity in 
Eq . (21) for the solution, while we find the cloud amount in Eq . (23). Let 
us tabulate the number distribution of positive and negative values of Apq 
for the solution x*l derived from 200 samples at each cloud amount and 
optical depth combination (Table 1) . We then further convert the number 
distribution to the percentage distribution by dividing the numbers of 
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Table 1. &oq Distribution - Number and oercentaae distribution for the oositive and 
neaative values of ioa for cases x»l in 200 samoles at eaor cloud amount and 
ootica! deoth combination 


H 

0 . 

(+) 

5 

(-) 

1 

(+) 

.0 

(-) 

1 . 

/ 

(+) 

5 

(-) 

2 . 

(+) 

0 

(-) 

2 . 

(+) 

5 

(-) 

3 . 

(+) 

,0 

(-) 

0.1 

8 

52 

14 

51 

19 

o 

52 

17 

68 

13 

67 

12 

78 

( Z) 

13 

87 

22 

78 

27 

73 

20 

80 

16 

84 

7 

93 

0.2 

19 

35 

18 

36 

15 

42 

16 

54 

8 

66 

6 

78 

(2) 

35 

65 

33 

67 

26 

74 

23 

77 

11 

89 

7 

93 

0.3 

32 

24 

25 

19 

6 

22 

4 

41 

0 

48 

0 

61 

(%) 

57 

43 

57 

43 

21 

79 

9 

91 

0 

100 

0 

100 

0.4 

36 

13 

13 

2 

4 

11 

2 

27 

0 

33 

0 

54 

( Z ) 

74 

26 

87 

13 

27 

73 
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either positive or negative value of Apq by their sum and multiplying 100 
as percentage. Results of the percentage distribution for the positive and 
negative values of Apq are plotted in Fig. 6. 

In Fig. 6, the percentage distributions of positive values of Apq are 
denoted by the numbers without parenthesis, and those of negative values 
are marked by the numbers with parentheses. Since the total percentage 
must be 100%, all the contours are shared by both positive and negative 
percentages. The areas without coverage of contours are the sections which 
pose no controversy of obtaining x“L in cloudy atmospheres . It shows in 
Fig. 6 that the higher percentage of positive values generally distributes 
in the areas of thinner cloud (or lower emissivity); while the negative 
values distribute more in the areas of thicker cloud. Furthermore, the 
gradient of positive value distribution is parallel to the abcissa (cloud 
amount) , but the gradient of negative value distribution is oriented along 
with the ordinate (cloud optical depth). Therefore, it is desired to find 
the cloud amount directly if Apq is positive, and find the cloud emissivity 
if Apq is negative. 

b. Application of the cloud retrieval program to HIRS/2 channels. 

The NOAA 6 satellite was launched on 27 June 1979. It is a third 
generation operational, sun-synchronous, near polar orbiting satellite with 
a height of about 810 km and having morning descending or afternoon 
ascending. Successive orbit3 cross the equator with 98.74 degrees of 
inclination and the orbital period is about 101.26 min. 






Page 24 


The NOAA 6 HIRS/2 instrument is an adaptation of the HIRS/1 instrument 
designed for and flown on the Nimbus 6 satellite. This instrument, being 
built by the Aerospace/Optical Division of ITT, measures incident radiation 
in 20 spectral regions of the infrared spectrum, including both longwave 
(15 pm) and shortwave (4.3 um) regions. The instrument cross-track scans 
+49.5 (+1120 km) to satellite subtrack. There are 56 scan steps per scan 
line with a resolution of 17.4 km near nadir and 29.9 km at the extremes of 
the scan. 

The data that are routinely processed at NOAA/NESDIS are available on 
nine-track, 1600 bpi tapes. The data of each orbit were packed in such a 
way that 20 channels of data were stored separately in 20 files. These 
data have been rearranged to contain co-located radiance values for all 20 
channels. The data set to be used in this analysis will be 24 April 1980, 
during which all channels of the HIRS instrument were operating properly. 

In numerical computations , the simulated atmosphere will be divided in 
such a manner that it coincides with the presssure levels used in the clear 
column radiance program (CCR) developed at NOAA/NESDIS. There are 40 
pressure levels for the OCR program. The program utilizes predetermined 
transmission profiles which can be empirically adjusted as a function of 
the known temperature and moisture profiles. The transmittance programs 
were kindly provided by L. McMillin. The weighting functions of the HIRS 
channels are depicted in Fig. 1. We will use three SIRS channels including 
two in the 15 pm CO^ band for cloud height retrieval and combining 15 um 
band wing channel and the 11 um window channel for cloud amount and 
emissivity calculations. The retrieved cloud parameters will be compared 
with available surface observations and satellite pictures to assess the 


Page 25 


reliability of the proposed cloud detection technique utilizing the 
existing data. In the course of this investigation, a concept of mapping 
various clouds over the globe may be developed and proposed for future 
satellite experimentation. 

c. Results of numerical calculation. 

Results of the sequential calculations for the cloud top height, cloud 
amount and emissivi ty are plotted in Figs. 7, 8 and 9. The figures show 
the mean value (bias) and standard deviation (S.D.) of 200 samples for each 
cloud amount and optical depth combination. The measurement and Instrument 
errors, which are created randomly as described in the last section, are 
imposed onto the exact "measured" radiances for the theoretical 
calculations . 

First, we assume the cloud height is at 9.37 km, which is 
approximately equivalent to 300 mb for the midlatitude spring 
climatological profile. By giving various combinations of cloud amount and 
optical depth, we shall examine the bias and stability of the alogrithm. 

The probability of this alogrithm miscalculating a cloudy atmosphere as a 
clear atmosphere or vise versa is shown in Fig. 7a. In addition to the 
criterion as described earlier to differentiate clear and cloudy 
atmospheres if x«l, we set another criterion such that if either AG^ or AG^ 
in Eq . (9) is equal to or less than 0 , the atmosphere is assumed to be 

clear. We found that the probability of mistreating a clear atmosphere as 

%• 

a cloudy atmosphere is generally less than 2%. However, for cloudy 
atmospheres with cloud amount equal to 0.1, we have about 38% of chance 
miscalculating them as clear atmospheres for t=* 0.5 or 7% for t= 3.0. The 
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cloud-free percentage clearly decreases with increase of cloud optical 
thickness even if the optical depth is greater than 1.5. This, of course, 
is the result of the treatments reevaluating the solution x»l as described 
in the last section. When the fractional cloud cover increases within 
IFOV, the cloud-free probability decreases rapidly, especially for 
optically thick clouds. 

For the cases we have determined as cloudy atmospheres , we then 
calculate the mean cloud height by averaging the solutions of those cloudy 
cases. Results show that the mean cloud heights are generally 
overestimated (Fig. 7b). This is particularly evident for optically thin 
clouds and/or small fractional clouds within IFOV. The S.D. shows that the 
stability of the cloud height retrieval is less satisfactory for optically 
thin clouds and the small fractional clouds. The S.D. decreases 
significantly with the increases of cloud amount and optical depth. 

The mean cloud amounts are estimated in such a manner that the results 
of the computations for the entire 200 samples are used for computing the 
mean values and S.D. (Fig. 7c). This means that for cloudy atmospheres 
which are miscalculated as clear atmospheres in the algorithm, results of 
A=0 will be counted in the estimations of mean cloud amount and the S.D. 

The outcome is as expected that the mean cloud amounts for optically thin 
clouds are underestimated. This is because the larger percentage of 
optically thinner clouds are mistreated as cloud-free in the alogrithm. 
However, the cloud amounts are generally overestimated if T>1. The outcome 
of the S.D. also shows that the algorithm is less stable when a optically 
thin cloud is present in the atmosphere. The underestimation of the cloud 
amount is also found for the case which is assumed to be nearly overcast . 
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This Is result of the constraint we impose to limit the solution of A to be 
between 0 and 1. Overall, the results are impressive, because the largest 
bias of the mean calculated cloud amounts is no worse than -0.18 and the 
S.D. is no greater than 0-25 for the cases of 0.5 <t<3.0. 

In Fig. 7d, the mean values of the calculated cloud emisslvltles and 
their S.D. are shown. The samples which are counted in for estimating the 
mean values and S.D. are those determined as cloudy atmospheres in the 
algorithm. The trend of the calculated cloud emissivity is, in general, 
opposite to that of calculated cloud amount. The bias of the calculated 
cloud emissivity is positive for optically thin clouds, but is negative if 
t> 1. The underestimation of the cloud emissivity for t>2 is obvious in the 
small fractional cloud cover areas. This may be because the gradient of 
percentage distribution for negative Apq is not distinct when the cloud 
amount is small. If we review Fig. 3, we find that when the cloud amount 
is small, the value of Apq is indifferent with respect to the variations in 
opical depths. 

Next, we test the algorithm for the different cloud heights. We find 
that for the cloud height assumed to be at 500 mb (-5.72 km), the 
probability of mistreating a cloudy atmosphere as clear atmosphere 
decreases compared to that of cloud height at 300 mb. The probability of 
the cloud-free percentage for cloud amount equal to 0.1 is 262 when t= 0.5 
and less than 12 when 1=3.0 as shown in Fig. 8a. In contrast to the 
conclusions by Wielicki and Coakley (1981), we find that the solution of 
cloud height retrieval for the lower cloud (at 500 mb) is not worse than 
the result of higher cloud (300 mb) . The accuracy of the cloud height 
retrieval is dependent on the degree selected for a polynomial describing 
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the relation between H function and Z £ in Eq. (9). By using lower degree 
such as 4 for the polynomial, the cloud height retrieval for the higher 
cloud may have a better solution than that for the lower cloud. This is 
because the polynomial of lover degree is better describing the relation of 
H vs. Z t for higher cloud. However, the polynomial of higher degree is 
more suitable for retrieving the lower cloud. The results of the cloud 
height retrieval based on the polynomial of different degrees for different 
cloud heights are shown in Table 2. Without imposing measurement noises in 
the experiment, the inherent bias of the polynomial is shown in the Table. 
In this study, the polynomial of degree equal to 5 is selected for the 
calculation. 

Table 2. Cloud Height Retrieval Bias : Sias of the cloud heioht retrieval based on the 

polynomial of different degrees for different cloud heiahts. 


N^egree 

Z t (km) 


5 

6 

9,370 

-0.089 

-0.173 

-0.239 

5.720 

0.120 

-0.069 

-0.104 

3.078 

-0.315 

-0.005 

-0.008 


In Fig. 8.c, the biases of the mean cloud amounts and their S.D. for the 
cases of cloud height at 500 mb are shown. In general, the patterns of the 
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bias and S.D. are quite similar to those at 300 mb height. Since the 

better solutions of cloud height retrieval result in better estimations of 

p and q in Eq . (12), and R and R in Eqs • (20) and (22), respectively, 

m n 

thus the solution of the cloud amount are also improved. The improvements 
are more evident in the regions of thinner clouds. Finally, the mean value 
and S.D. of calculated cloud emissivity of cloud height at 500 mb have the 
same trend as those of cloud height at 300 mb, also with slight improvement 
of the solutions in thinner clouds region. 

For the cloud height at 700 mb (~3.08 km), the results of cloud-free 
percentage, cloud height, amount and emissivity are shown in Figs. 9a, b, c 
and d, respectively. The general patterns of the solutions are quite 
similar to those of cloud heights at 300 and 500 mb. Overall, the 
solutions for lower clouds are as satisfactory as those for higher clouds . 
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Fid. 9a Cloud Amount versus Optical Oeoth (restrictive criteria imposed) for Cloud Ht at 
700 mo 
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4. Applications to HISS/2 data of NOAA 6 and TISOS-N Satellites 

On the basis of the previous depiction and analysis it is evident that 
a proper combination of IS sounder data may be successfully used to recover 
the cloud parameters.lt appears desirable and Important to apply the 
developed retrieval scheme to real satellite data to investigate its 
applicability. Total number of 28 cases with different synoptic situations 
were chosen to perform the retrieval exercises. Results are illustrated 
and will be discussed in this section. 

Since NOAA 6 and TISOS-N pass the Great Plains at 1400 and 1000, 
respectively, the mean value of colocated radiances measured by both 
satellites may closely represent the observed data at 1200, at which time 
Che SA0B data are available. The RA0B profiles are employed as input data 
for the calculation of clear atmospheric radiances, and the surface reports 
are used for the comparsion with the computed results. 

There were cases when the RAOB stations did not coincide exactly with 
HISS footprints. Because of the horizontal variability of the HISS data, a 
distance-weighted method has been devised to obtain an average measured 
radiance from the available HISS data around the RAOB station. This 
procedure is especially important for wing channels which are affected 
significantly by the surface condition. All of the radiosonde observations 
used in this study were available up to at least 100 mb. 

In this paper, we will present cases on 24 April 1980 along with the 
GOES-East satellite IS picture (Fig. 10) and the corresponding synoptic 
observations (Fig. 11). No GOES visible picture will be used, because the 
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Fiq. 10 GOES Ifi Imagery for April 2A, 1980 - imagery is for 1200 GMT with the frontal 
position shown. 



Fiq. 11. Computed Cloud Height Analysis Analysis from HIRS/2 data is shown alono with 
Svnootic reports for 1200 GMT 2A April 1980 
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Great Plains at 1200 was still in dark. On 24 April, there was a low 
pressure in the eastern Great Plains. The cold front stretched out from 
the low center to northern Texas , and stationary front extended into the 
Atlantic Ocean. Analyses of surface reports and GOES image Indicate that 
cirrus clouds were present in the stationary front area, while multi-layer 
clouds and/or precipitation occurred in the vicinity of the cold front 
area. 3ased on the availabity of RAOB profiles, 28 stations were selected 
for the retrieval exercises . The surface reports of the chosen stations 
are recorded as shown in Fig. 11, and the results derived from the 
numerical computations are also plotted. 

Results of the cloud height, amount and emissivity retrieval indicate 
that a strong convective cloud system may exist in large portion of 
Oklahoma and Kansas states. The high cloud top height, which were found to 
be as high as 10 km (Fig. 11), and the high emissivity in this area 
(Fig. 13 ) suggests that the system is either composed of thick clouds or 
multi-layer clouds. The cloud amount in this region was close to 
overcast as shown in Fig. 12. 

This convective cloud band is obviously illustrated in the GOES IR 
picture (Fig. 10), in which the positions of frontal systems are marked. 
While in the southeast of this cloudy area (eastern Texas and western 
Lousiana) , extensive low clouds were found according to our retrievals . 
Although these low clouds were not clearly seen from the IR picture, the 
surface reports confirmed the existence of the clouds. The large values of 
the cloud amount and emissivity, and the low cloud top height indicate that 
this might be a less convective area compared to the post-frontal area. 

The low clouds might have reached the ground, and looked more like fog from 
some of the surface reports. To the north of the stationary front, ex- 
tensive cirrus clouds were found. The cloud top heights were at least 
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9 km, and Che emissivities were about 0.6. The cirrus clouds in West 
Virginia slight have some scattered low clouds underneath because large 
values of cloud emissivity (>0.7) were found. There were also large areas 
of low clouds existing in the southeast of the Great Lakes covering 
Michigan and northern parts of Illinois and Indiana. To the south of the 
stationary front, there was largely clear atmosphere with only small 
fraction of high clouds detected. Those high clouds could be thin cirrus 
because of small emissivities as the calculated results show. 
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5. Conclusion 


In Chls study, a numerical method for computing cloud top height, 
amount and emissivity has been developed based on the parameterized 
infrared radiative transfer equation for cloudy atmospheres. Theoretical 
studies were carried out based on synthetic atmosphere containing various 
thickness of fractional cloud cover within a EFOV. Upon imposing 
instrument noises to the upwelling radiances the solutions of the cloud 
parameters under different cloud conditions are obtained. It is found, that 
the measurement error effect is more profound for the thinner and/or 
smaller fractional cloud case in the numerical calculation. But, in 
general, resulting cloud parameters retrieved from the present program are 
found to reasonably accurate by comparision with theoretically assumed 
solutions. Numerical experiments were also performed utilizing NOAA 6 and 
TIB.OS-N HIRS/2 data in which a mesoscale cloud system is analyzed. 

Numerical results show that the cloud areas were described distinctively in 
terms of the cloud top height, amount and emissivity derived from the 
program- 

Moreover, the computed cloud height and amount are generally in good 
comparison with those obtained from synoptic reports and satellite picture. 
However, further verification of the satellite cloud sounding techniques 
requires carefully designed field experiments in which reliable cloud 
properties, such as cloud top height, phase, thickness and emissivity, may 
be obtained from aircraft observations under the satellite pass. 

In application of the current technique to achive cloud information in 
large scale study (e.g., ISCCP), the satellite radiances have to be 
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arranged In Che matrix manner. The entire area within the scale needs to 
be divided into smaller grid boxes (such as 2. 5x2. 5°), in which the 
radiance histogram can be constructed. The clear column radiance within 
the grid box can be derived based on the histogram. This clear column 
radiance may then be used in this cloud retrieval program to replace the 
clear column radiance computed based on the radiosonde profiles in this 
study. Such arrangement may significantly save computing time on analyzing 
large amounts of satellite data. Further study concerning more efficient 
cloud retrieval technique in the large scale is under investigation and 
will be reported in future. 


- H'*,. iUi 
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